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Abstract 

In this note, we announce a general result resolving the long-standing question of nonlinear 
modulational stability, or stability with respect to localized perturbations, of periodic traveling- 
wave solutions of the generalized Kuramoto-Sivashinski equation, establishing that spectral 
modulational stability, defined in the standard way, implies nonlinear modulational stability 
with sharp rates of decay. The approach extends readily to other second- and higher-order 
parabolic equations, for example, the Cahn Hilliard equation or more general thin film models. 



1 Introduction 

In this note, we describe recent results obtained using techniques developed in [JZ) IJZN|, IBJNRZ2| 
on linear and nonlinear stability of periodic traveling-wave solutions of the generalized Kuramoto- 
Sivashinsky (gKS) equation 

(1.1) ut + jd^u + ed^u + Sdlu + d^{u^ /2) = 0, S > 0, 

a canonical model for pattern formation in one spatial dimension that has been used to describe, 
variously, plasma instabilities, flame front propagation, turbulence in reaction-diffusion systems, and 
thin film flow down an incline (SH |KTl EHl iDSSllPSU] . 

More generally, we consider (taking without loss of generality 7 = 1) an equation of the form 

(1.2) ut + dlu + edlu -t- 5dlu + d^f{u) = 0, 

f G C^, S not necessarily positive. Our methods apply also, with slight modifications to accomodate 
quasilinear form (see |JZN| ) to the Cahn Hilliard equation and other fourth-order models for thin 
film flow as discussed for example in jBMSZj . Indeed, the argument, and results, extend to arbitrary 
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2r-order parabolic systems, so is essentially completely general for the diffusive case. As shown in 
[JZj . they can apply also to mixed-order and relaxation type systems in some cases as well. 



It has been known since 1976, almost since the introduction of the model (1.1) in 1975 |KT1 ISl) . 
that there exist spectrally stable bands of solutions in parameter space; see for example the numerical 
studies in |CKTR[ IFSTj . Moreover, numerical time-evolution experiments described for example in 
[CD| suggest that these waves are nonlinearly stable as well, serving as attractors in the chaotic 
dynamics of (gKS). However, up to now this conjecture had not been rigorously verified. 

Here, we announce the result, resolving this open question, that spectral modulational stability^ 
defined in the standard sense of the modulational stability literaturej^ implies linear and nonlinear 
modulational stability. Our analysis gives at the same time new understanding even at the for- 
mal level of Whitham averaged equations. Further details, along with numerical investigations of 
existence and spectral stability, will be given in [BJNRZ3) . 

2 The traveling-wave equation 



Substituting u — u{x — ct) into ( |1.2[ ), we find that the traveling-wave equation is —cu' + u"" -\- eu" 
6u" + /(u)' = 0, or, integrating once in x, 

(2.1) -CU + u'" + eu" + Su' + f{u) = q, 

where g G M is a constant of integration. Written as a first-order system in {u,u' ,u"), this is 
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(2.2) 
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It follows that periodic solutions of ( |1.2[ ) correspond to values {X, c, q, b) S M^, where X, c, and q 
denote period, speed, and constant of integration, and b = (&i, 62, ^3) denotes the values of (m, m', u") 
at a; = 0, such that the values of (m, u', u") at a; = X of the solution of \2.2\ are equal to the initial 
values (&i,52,63). 

Following [JZ], we assume: 

(HI) / e C^+\ if > 4. 

(H2) The map H -.W^ ^M? taking {X, c, q, b) ^ {u, u', u"){X, c, b;X)-b is full rank at (X, c, 6), 
where {u,u' ,u"){-] •) is the solution operator of ( |2.2| ). 

By the Implicit Function Theorem, conditions (H1)-(H2) imply that the set of periodic solutions 
in the vicinity of U form a smooth 3-dimensional manifold (counting translations) 

(2.3) {U'^{x - a - c(/3)t)}, with aeR, f3 eR^. 

3 Bloch decomposition and spectral stability conditions 

In co-moving coordinates, the linearized equation about u reads 

(3.1) vt = Lv := ((c - a)v)^ - v^xxx - ev^xx - Svxx, a := df{u), 

and the eigenvalue equation as Lv := —Vxxxx ~ ^i^xxx ~ ^v^x + ((c — '^'v) ^ = ^"^^ Following [U], we 
define the one-parameter family of Bloch operators 

(3.2) :=e-*«"Le*«", $ G [-tt, tt) 

1 In particular, the sense verified in ICKTRI IFST] 



4 WHITHAM AVERAGED EQUATIONS 



3 



operating on the class of periodic functions on [0,X]; the (i^) spectrum of L is equal to the 
union of the spectra of all with ^ real with associated eigenf unctions w{x, ^, A) := e^^^q{x, ^, A), 
where q, periodic, is an eigenfunction of i^. By standard considerations, the spectra of consist 
of the union of countably many continuous surfaces Aj(^); see, e.g., jG]. 

Without loss of generality taking X = I, recall now the Block representation 

(3.3) uix)^[^) J\^i-^ui^,x)d^ 

of an function u, where u{£,,x) := J^k e'^^^^^u{^ + 27rfc) are periodic functions of period X = 1, 
u{-) denoting with slight abuse of notation the Fourier transform of u in x. By Parseval's identity, 
the Bloch transform u{x) — > u{£^,x) is an isometry in L^: ||u||^2(^-) = Jl^ Jq dx — 

\\u\\l^(^;L^{x))j where L'^{x) is taken on [0,1] and L^(^) on [— 7r,7r]. More generally, for q < 2 < p, 
i + i = 1, there holds the generalized Hausdorff- Young's inequality [JZj 

(3.4) I1'"IUp(x) < \\u\\Llii:LP{x))- 

The Bloch transform diagonalizes the periodic-coefficient operator L, yielding the inverse Bloch 
transform representation 

(3.5) 6^*^.0 = (^) f e'^--e^^'uoi^, x)d^ . 



Following [JZj . we assume along with (H1)-(H2) the strong spectral stability conditions: 
(Dl) am) C {ReA < 0} for C 7^ 0. 

(D2) Recr(L5) < -6l|^|2, 6* > 0, for ^ e M and |C| sufficiently smaU. 
(D3) A = is an eigenvalue of Lq of multiplicity 2j^] 

Assumptions (H1)-(H2) and (D1)-(D3) imply |JZllBJNRZ3j that there exist 2 smooth eigenvalues 
(3.6) A,(0 = -*a,e + o(|ei), J = 1,2 

of bifurcating from A = at ^ = 0. Following jJZj . we make the further nondegeneracy hypotheses: 
(H3) The coefficients aj in (3.6) are distinct. 

(H4) The eigenvalue of Lq is nonsemisimple, i.e., dimkerLg — 1- 



4 Spectral stability and the Whitham averaged equations 



As noted in [Sell JZj . coefficients Oj in (3.6) are characteristics of the 2x2 Whitham averaged system 

Mt + = 0, 
(4-1) , , 

LOt + {C^)x = 

formally governing large-time (~ small frequency) behavior, evaluated at the values c, w of the 
background wave u, where M(c, w) is the mean of u over one period of the periodic wave with speed 
c and frequency w = 1/A', and F{c,lu) is the mean of f{u). Here, uj ^ ipx, c ^ —ipt/i'x, where 'ip 
denotes phase in the modulation approximation 

(4.2) u{x,t) u{i;{x,t)). 



^ The zero eigenspace of Lq, corresponding to variations along the 3-dimensionaI manifold of periodic solutions in 
directions for which period does not change IJZI . is at least 2-dimensional by (H2). 
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In the context of ( 1.1 1, thanks to the Galillean invariance x ^ x — ct, u ^ u + c, (4.1 1 reduces to 
3^ ct + {H{uj) - to(w)c)^ = 0, 

LOt + icuj)x = 0, 

where m(a;) denotes the mean over one period of u for a zero-speed wave of frequency w, and H{uj) 
the mean of and in the classical situation e — considered in [FSTj . to Ct + {H{uj))x — 0, 

ujt + {cio)x — 0, which hnearized about background values c = 0, uj = wq, yields a wave equation 

(4.4) V« + iOQdHiLUo)il^xx = 

so long as dH(ujo) < 0. Indeed, by odd symmetry, we may conclude in this case that the second-order 
corrections bj in the further expansion 

(4.5) X^{^)^^a,^^b,e■■■ 



of (3.6) are equal, hence Aj(^) agree to second order with the dispersion relations of the viscoelastic 



wave equation 

(4.6) iptt + ^^odH{ujo)ip^^ ^ d{oJo)il;txx, d =261 = 262- 

This recovers the formal prediction of "viscoelastic behavior" of modulated waves carried out in 



[FST| and elsewhere, or "bouncing" behavior of individual periodic cells. Put more concretely, (4.6) 
predicts that the maxima of a perturbed periodic solution should behave approximately like point 
masses connected by viscoelastic springs. However, we emphasize that such qualitative behavior- 
in particular, the fact that the modulation equation is of second order- does not derive only from 
Gallilean or other invariance of the underlying equations, as might be suggested by early literature 
on the subject, but rather from the more general structure of conservative (i.e., divergence) form 
[SellJZ| F] Indeed, for any choice of /, XjiO may be seen to agree to second order with the dispersion 



relation for an appropriate diffusive correction of (4.1 ), a generalized viscoelastic wave equation. See 



[NRll INR2) for further discussion of Whitham averaged equations and their derivation. 



5 Linear estimates 

The main difhculty in obtaining linear estimates is that, by (D3) and (114), the zero eigenspace 
of Lq has an associated 2x2 Jordan block. This means that e^"* is not only neutral but grows 
as 0{t). Viewed from the Bloch perspective, it means that the eigenprojections of blow up as 
^ — > 0. Performing a careful spectral perturbation analysis, separating out the singular part of the 



eigendecomposition of e « in (3.5|, and applying (3.4), as in Lemma 2.1, Prop. 3.3, and Prop. 3.4 



of |JZ] . we obtain the following detailed description of linearized behavior. 

Proposition 5.1. Under assumptions (H1)-(H4) and (D1)-(D3), the Green function G{x,t;y) of 



(3.1) decomposes as 



(5.1) G{x, t; y) = u'{x)e{x, t; y) + G{x, t; y), 



^As discussed further in Z,, conservation of mass lies outside the usual Noetherian formulation. 
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where, for some C > and all t > 0, l<q<2<p<oo and 1 < r < A, 
G{;t-y)f{y)dy 



(5.2) 
(5.3) 
(5.4) 



d;Gi;t;y)f{y)dy 



dtG(-,t-y)f{y)dy 



< Crii----p) (1 + <)-K|-^-^) ||/||^,nL^ 

1/1 1\ 7- 1/2 11^1|T- 

1/1 1\ ^ l/'2 1 l^il 

e{x,t;y) = forO <t <1, and for all t > Q , l<q<2<p<oo,0<j,l,j + l<K, and 1 < r < A, 
did\e{;t-y)f{y)dy 



LP( 



LP( 



(5.5) 



, 1 (' 1 1\ {j + k) 

<C{l + t) p) 



did\dle{;t;y)f{y)dy 



LP{ 



LP( 



1 1(1 i\ (j+fc) 



./IU''(R) 
ll/IU'fR 



Moreover, for some constants pj, and Oj and bj as in (4.5), 



(5.6) 



ei-,t]y) - ^Pjerrfn 




Defining ip := —e and noting that u{x)+ij^{x, t)u' (x) ~ u(x+Tj){x, t)), we see from (5.1 H( 5.6 1 that 
hnearized behavior indeed agrees to lowest order with modulation by a phase function ip satisfying 



a generalized viscoelastic wave equation obtained by diffusive correction of (4.1), consisting of a 
first-order hyperbolic-parabolic system in tpx, ipt. 



This observation not only generalizes the second-order scalar description (4.6) obtained in the 



special case e = |FST) . but is in some sense more correct. For, note that the second-order scalar 
description can be a bit misleading as regards the assumption of initial data. In particular, from 



the description (4.6), one might be tempted to conclude that the linear response to a compactly 



supported perturbation would consist of the D'Alembertian picture of two approximately compactly 
supported wave forms in ip moving in opposite directions, diffusing slowly at Gaussian rate. Yet, 
the explicit bound (5.6) shows that this is rather a description of the derivative il)x\ 



Indeed, as described further in Section 1.2, |JZN] . it is the variables V's, i^t that are primary, rather 
than ?A,'0t as suggested by (4.6 1, and it is these variables that are related to the initial perturbation 



{u — u)\t=o- Thus, our analysis gives not only technical verification of existing observations, hut also 
new intuition regarding the nature of modulated behavior. 

6 Nonlinear stability 



Using the linear bounds of Prop. |5.1| together with nonlinear cancellation estimates as in |JZ1 IJZN] . 
we obtain, finally, our main result describing nonlinear behavior under localized perturbations. 

Theorem 6.1. Assuming (H1)~(H4-) and (D1)~(D3), let u = (f,u) be a traveling-wave solution of 



pr2| ). Then, for some C > and e W^'°°{x,t), 

\\u-u{--^- ct)\\Lp{t) < C{1 + - u\\LinH'<\t=o, 

(6.1) \\u-u{-~i^~ ct)\\HK (t) < C{1 + t)--^\\u - u\\LinH'< |t=o, 

' '"^"u - "IlLinff'f |t=o, 



|(Vt,V'x)||w-+i.p < C(l -f i)-5(i-i/f)llf' - 



7 APPLICATION TO KURAMOTO-SIVASHINSKY 



and 
(6.2) 



\\u-u{- - ct)\\L'^{t), 11-0(01^°° < - M||i,in//A'|t^o 



for all t > 0, p > 2, for solutions u of (1.2) with \\u~ u\\]^ir]H'^'\t=o sufficiently small. In particular, 
u is nonlinearly bounded D — ^ L°° stable. 



Similarly as in the discussion of linear behavior, we note that Theorem |6.1| asserts asymptotic 
n — !• LP convergence of u toward the modulated wave u{x — tp{x,t)), but only bounded 
L^ n — )■ L°° stability about u(x), a quite different picture from that suggested at first sight by 



(4.6) 



7 Application to Kuramoto— Sivashinsky 

We conclude our discussion by displaying some representative traveling wave orbits and their as- 
sociated spectrum, for the case e = 0.2, computed respectively using MATLAB and the spectral 
Galerkin package SpectrUW. These indicate, similarly as in the e = case studied in [CKTRi FFSTI . 
the existence of a band of spectrally stable periodic traveling waves. For related studies, and an 
animation of the spectral evolution, see Gallery, [http: / /www. math. i ndiana.edu/ For more detailed 
numerical verification using the periodic Evans function |G ) . see |BJNRZ3] . 






Figure 1: Evolution of spectra as frequency w = 1/X varies with c = held constant. Here e = 0.2, 
(5 = 1,7 = 1. Running left to right in the top row, we see the evolution of the periodic orbits in 
the three-dimensional phase portrait starting from Hopf bifurcation and moving through the stable 
band of periodic waves, with frames directly below depicting the spectrum of the corresponding 
linearized operator about the wave. Orbits were computed with MATLAB; spectra were computed 
with the SpectrUW package developed at University of Washington |CDKK| . 
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